Phase diagram of cold atoms in the bilayer hexagonal optical lattices 
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We propose a scheme to realize the supersohd of cold atoms in the bilayer hexagonal optical 
lattices. We show that the bilayer optical lattices can be realized by coupleing two monolayer 
optical lattices with two vertical standing laser waves. Using the quantum Monte Carlo method, we 
calculate the ground state phase diagrams which contain an antiferromagnet, a solid, a superfluid, a 
fully polarized and especially, a supersolid. A supersolid phase emerge when the lattice anisotropy 
or the ratio of intra- to inter-layer tunneling increase to a critical value and it disappear when the 
anisotropy or the tunneling ratio are big enough. We give a experimental protocol to observe those 
novel phenomenon in the real experiments. 
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Introduction. — In recent years, ultracold atoms in op- 
tical lattices have been widely used to simulate many- 
body phenomena in a highly controllable environment 
[1-6]. By designing configurations of the atomic sys- 
tem, one can simulate effective theories of the forefront 
of condensed-matter physics. Recent researches show 
that graphene have fascinating effects and exhibit par- 
ticularly rich quantum phases [7-12]. Compared with 
the graphene, the systems with cold atoms in an optical 
lattices are much more controllable. Very recently, the 
realization of the tunable spin-dependent hexagonal lat- 
tices indicate that it can investigate the quantum phases 
of many systems which has a hexagonal geometry [13, 14]. 

One of the goals of studying the cold atoms in optical 
lattice is to search for the novel states. Supersolid phase 
(SS) is an exotic state, where both the diagonal and off- 
diagonal long range order coexist. Although this novel 
state is still not find in experiments [15-19], its pres- 
ence has been confirmed theoretically in lattice model, 
including bosonic systems [20-28] and spin systems [29- 
33]. As a powerful tool to tailor the quantum phases, it 
is reasonable to ask whether it can realize the SS phase 
in optical lattices ? Compare to the monolayer material, 
the bilayer systems show fascinating different properties 
[34, 35]. Can we design a tunable bilayer hexagonal op- 
tical lattices base on the monolayer lattices ? 

In this Letter, we propose a scheme to realize the su- 
persolid in bilayer hexagonal optical lattices. As shown 
in Fig. 1, the bilayer hexagonal optical lattices can be 
formed by coupling two monolayer hexagonal lattices 
with two vertical standing waves, where one standing 
wave has the twice period of the other. The monolayer 
lattices can be set up by intersection of three laser waves 
at an angle of 120°. Using stochastic series expansion 
(SSE) quantum Monte Carlo (QMC) method [36] we cal- 
culate the ground state phase diagram of cold atoms in 
this bilayer optical lattices. Results show that the SS 
can be realized by adjusting the lattice anisotropy or the 
intra- and inter-layer tunneling and this bilayer optical 



lattices could be an ideal instrument to tailor the su- 
persolid state. We give a protocol of the observation of 
supersolid in real experiments at last. 




FIG. 1. (Color online) The upper plot show the schematic dia- 
gram of the formation of the bilayer hexagonal optical lattice. 
The yellow and cyan planes denote the two monolayers of the 
lattice which were formed by interacting of three standing 
laser waves. The two big arrows show the two vertical stand- 
ing waves used to couple the two monolayers. Red sphere 
(green sphere) denote sublattice A (B). The thin line denote 
the intra-plane exchange coupling t„ and the blue doublear- 
row denote the inter-plane exchange coupling t'^. The lower 
left plot show the top view of monolayer optical lattice. In the 
lower right plot, the yellow uparrow and magenta downarrow 
denote the two relevant internal states of one atom. U^, U_i 
correspond to the strength of the on site interaction energies 
for the cold atoms that have two same states ^ t and J, 4, f/^ j. 
is the interaction strength for the cold atoms that have two 
different states t4- 
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FIG. 2. (Color online) Ground state phase diagram of — 
h for Eq. 1 with J = 0.3. The symbol is the Monte Carlo 
calculations and solid (dashed) lines denote the boundaries of 
the first (second) order phase transition. The phase diagram 
contain an antiferromagnet (AF), a superfiuid (SF), a solid, a 
supersolid (SS) and fully polarized (FP). 
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FIG. 3. (Color online) The evolution of the order parameters 
pb, Sa and magnetization with h increasing for different 
Jz. In the calculation, L = 8, 10, 12, 16 and J = 0.3 are used. 



The bilayer hexagonal optical lattices. — Considering 
cold bosonic atoms (^^Rb) in an optical lattice, we as- 
sume that each atom has two relevant internal states 
(F = 1, mp = ±1); which are denoted by the two 
spin index a =t, i and the atoms with a =t, i are 
trapped with independent standing wave laser beams 
through polarization selection. The three dimensional 
hexagonal lattice can be set up by intersection of three 
coplanar laser beams under an angle of 120° between 
each other in the x-y plane and two intersecting waves 
along the z direction. The total potential of the lattice 
is U{x,y,z) = Y,i=i,2,3^cr^^'"^l^Lixcos9i + ys'mdi) + 
TT/2] + VziaCos'^{kLz) + Vz2aCOs'^{2kLz), wherc 9i = tt/S, 
6'2 = 27r/3, 6*3 = and — 27t/X is the optical wave 
vector with A — 830 nm. Va is the barrier height of a 
single standing wave laser field in the x-y plane [37, 38], 
Vzia and Vz2a are the barrier heights along z axis [3]. To 
form decoupled bilayer lattices, 4Vz2a > Vzia is required 
and the ratio of inter- to intra-bilayer potential heights 
should be so large that the tunneling between two neigh- 
bor bilayers can be prevented. 

Around the minima of the potential harmonic approx- 
imation has been taken and the spin-dependent bilayer 
optical lattice can be described by a Bose- Hubbard model 

H = - (taai^^aaja + H.c.) - Y{t'^a\^^a2^a + H.c.) 

a(i.j) iiJ 

aia ai 

~ l^'^inatt - natl) , (1) 
ai 

where a denotes the layer index 1,2, aj^ (ata) is the cre- 
ation (annihilation) operator of the bosonic atom at site i 



for spin cr, rii^ — aj^aia- is the occupation number on site 
i. runs over nearest neighbors. We will only focus 

on the regime of strong coupling, Ua, U^i ^ t<T, t'^, i.e. 
the mott insulator regime, where the atomic density of 
one atom per lattice site. 

The phase diagram of ground state. — Let J± = 

2t'/jUu, J = 2t^h/Utu Jz = -m^ + h^)/mh) + 

2Un[t^/{U^t^)+tJ{U^t^)] and h = 2ii-2{t\^ +t^-') IU^+ 

2{t\^ + h^)/Ui- In the following, we will take Jj_ = 1 
as energy unit and use the stochastic series expansion 
quantum Monte Carlo method to calculate the ground 
state phase diagram of Eq. 1. In our simulation, periodic 
boundary condition is imposed and the temperature is set 
inversely proportional to 2L with L — 8, 10, 12 and 16. 
The lowest temperature can reach to 0.03125, which is 
low enough to insure that we can investigate the ground 
state properties. 

In order to focus on the low energy physics, J ^ 1 is 
required (in this Letter 0.1 < J < 0.4). The four eigen- 
states of each two atoms coupled by the inter-layer tun- 
neling: sing let state (|s) = ^(| tl> - I It))) and triplet 

states (IM = ^(|n) + lit)),|ti) = ltt>,|t-i) = lU)), 
are conveniently used to describe the ground state of the 
system. When there exist a weak magnetic field, the 
triplets split into three and the lowest state (|ii)) keeps 
dropping with the field increasing. It is clear that the 
contribution of |to) and to the ground state is neg- 
ligible. When the magnetic field exceed a critical value, 
{\ti)) take the place of \s) to be the lowest energy level. 
A triplon bosonic operators fej = ( — l)*/\/2(s]^j — s^J 
can be used to transform \s) into (|ti)) [30]. It has the 
following operation: 6| |t_i)i = \s)i, bl\s)i = \ti)i and 
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FIG. 4. (Color online) Plots of (a) the number densities of the 
states: |s), \ti) and |to), (b) the SF and solid order parameter 
Pb and Sa and (c) the magnetization niz and the staggered 
magnetization ml changing with magnetic field h increasing 
for Jz = 3.3 with J = 0.3. 




FIG. 5. (Color online) (a) Finite size extrapolation of pt 
and Sa- The symbol is the QMC results and magnetic fields 
change from 2.1 to 2.2 with increments 0.02 from bottom to 
top for Pi, and from top to bottom for Sa- (b) The results of 
Pb and Sa for L=8, 10, 12 and 16. 



K \tih = 0. 

To character the supcrfluid state, we compute the 
condensate density pi, ~ ^/^'^J^iji^l^j) ^^le 
number of inter-layer bond, N ^ 2 x L^). For solid 
state, the usually static structure parameter S'(Q) is 
not suitable for the honeycomb lattice and we define 
the structure parameter Sa = 

where = +1 (—1) for i on sublatticc A (B), n.i rep- 
resent the number of (|ti)) on dimmer i. The antiferro- 
magnet is characterized by the staggered magnetization 
Ml = l/(2N){J2ieiiS(, - Sli))- The magnetization is 
calculated by = l/(2iV)(E,(5f, + ^f,))- 

The ground state phase diagram of H in the Jz — h 
plane with J = 0.3 is shown in Fig. 2. Order parame- 
ters Pb, Sa and vs the magnetic field h for different 
Jz are shown in Fig. 3. When Jz = 1, i.e., the isotropic 
case, the system is in the SF state, which is character- 
ized by > and 5*^ ~ 0. In this state, the U(l) 
symmetry of the system is broken by the Bose-Einstein 
condensates of \ti). When Jz > 1, without loss of gen- 
erality, we select Jz ~ 3.3, the evolution of the number 
densities and order parameters with the magnetic field 
h increasing is shown in Fig. 4. At lower magnetic field, 
the only nonvanishing order parameter show that the 
system is in the antiferromagnet state. This state is the 
traditional Ising order and the system has the staggered 
up and down arrangement in both intra- and inter-layer. 
When the magnetic field increased to a critical value, the 
staggered magnetization drops abruptly while the mag- 
netization has the inverse behavior, this denote the first 
order transition from AF to SF. Between 1.0 < Jz < 3.0, 
with magnetic field increasing, the system only experi- 



ence AF, SF and FP states. 

As the anisotropy Jz > 3.0, see Fig. 4. with the magnetic 
field increasing, the number density of |ti) jump abruptly 
to a finite value in the transition from AF to SF. Then, 
the triplet state keeps increasing to half filled and 
the system is in a solid state which is made of staggered 
singlet state and triplet state \ti). This state is charac- 
terized by Sa > and pb = and its obvious feature is 
the gapped plateau in the curve. After the platform, 
TT-ti keeps on increasing and the SF state reappear. In 
the intermediate region between SF and solid or between 
solid and SF, there appears a supersolid state which is 
characterized by both Sa > and pf, > in the thermo- 
dynamic limit. Finite size extrapolation results shown in 
Fig. 5 indicate that the SS is nonzero when L is extrap- 
olated to infinite. It is obviously that the transitions of 
SF-SS, solid-SS, SS-solid and SF-solid are second order. 
When Jz = 4.8, the SS state disappear, the transitions 
of AF-solid and solid-SF are first order. 

In the following, we want to consider the ground state 
phase diagram in the J — h plane at fixed Jz . The phase 
diagram with Jz = 3.3 is shown in Fig. 6. In the small 
J, the system only appear AF, SF, solid and FP. When 
J increase to nearly 0.15, two sharp regions of SS appear 
between SF and QS and they turn wider with J increas- 
ing and merger into one at J = 0.3, where QS phase 
disappear. Within 0.17 < J < 0.32 the SS state is stable 
in the thermodynamic limit and it vanishes and leaves 
SF phase only when J up to 0.35. 

Experimental protocol- — The experimental protocol of 
cold atoms in the bilayer hexagonal optical lattices can 
be taken as follows: *^Rb Bose-Einstein condensates were 
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FIG. 6. (Color online) Ground state phase diagram of J—h for 
Eq. 1 with the parameter Jz = 3.3. The symbol is the Monte 
Carlo calculations, solid (dashed) lines denote the boundaries 
of the first (second) order phase transition. The phase di- 
agram contain an antiferromagnet, a superfluid, a solid, a 
supersolid and fully polarized. 



created up to 10^ atoms in the state {F = 1, mF = — 1) by 
employing radio-frequency evaporation. Here F denotes 
the hyperfine states and mi? is the Zeeman substatcs of 
F. The cooled atoms were confined in a nearly isotropic 
crossed dipole trap with a trap frequency of w w 2irx. 90 
Hz. The preparation of different pure and mixed spin- 
states was performed with the aid of microwave sweep 
[3]. Then the optical lattice was ramped up in 80 ms. 

To form the three dimensional hexagonal lattices, three 
optical standing waves were aligned intersection under 
an angle of 120° in the x-y plane. Orthogonal to the x-y 
plane, two intersecting waves operated at the same wave- 
length were settled. One of the two waves had the twice 
period of the other. The laser beams were derived from 
a Ti:sapphire laser operated at a wavelength A = 830 nm 
(red detuned). The beams were linearly polarized within 
the intersecting plane [13, 14]. Using the optical fibers, 
all the beams were spatially filtered and guided to the 
experiments. 

To detect the magnetization and staggered mag- 
netization ml which in the range 0-0.5, a quantum po- 
larization spectroscopy method can be used: along the 
z axis a pulsed polarized light was sent to the atoms 
trapped in the lattice. After the light polarization cou- 
pled with the atomic spin, the atoms magnetization can 
be calculated by the light polarization rotates [39] . The 
superfluid density pi, has the maximum no more than 
0.25 could obtained from spin-spin correlation which can 
be detected by the noise correlation method [40]. Spin 
structure factor Sa has the maximum no more than 0.2 
can be calculated by the Bragg scattering [41, 42]. 

Conclusion. — In summary, we proposed a scheme to 
realize supersolid of cold atoms in a bilayer hexagonal op- 
tical lattices. We show that the bilayer hexagonal optical 



lattices can be set up by intersection of two intersecting 
waves orthogonal to the two monolayer hexagonal optical 
lattices. The intra- and inter-layer hopping and lattice 
anisotropy of the optical lattice can be designed by ad- 
justing of optical potentials and the additional zeeman 
field. The QMC calculation results show that this bilayer 
hexagonal optical lattices may be an practicable way to 
tailor supersolid. We hope that this optical lattices will 
be used to realize the supersolid in real experiments in 
the near future. 

This work was supported by the NKBRSFC un- 
der grants Nos. 2011CB921502, 2012CB821305, 
2009CB930701, 2010CB922904, NSFC under grants Nos. 
10934010, 11228409, 61227902 and NSFC-RGC under 
grants Nos. 11061160490 and 1386-N-HKU748/10 and 
the ISTCPC under grant No. 2012DFG11610. 



[1 
[2: 

[s: 

[4: 
[s: 

[6: 
[7: 

[8 

[9: 

[10: 
[11 

[12 
[13 

[14 

[15 
[16 
[17 
[18 
[19 



M. Greiner, O. Mandel, T. Esslinger, T. W. Hansch and 
I. Bloch, Nature 415, 39 (2002). 

N. Gemelke, X. B. Zhang, C. L. Hung and C. Chin, Na- 
ture 460, 995 (2009). 

S. Trotzky, P. Cheinet, S. Foiling, M. Feld, U. Schnor- 
rberger, A. M. Rey, A. Polkovnikov, E. A. Demler, M. D. 
Lukin and I. Bloch, Science 319, 295 (2008). 
C. Orzel, A. K. Tuchman, M. L. Fenselau, M. Yasuda, 
and M. A. Kasevich, Science 291, 2386 (2001). 
L. Hackermiiller, U. Schneider, M. Moreno-Cardoner, T. 
Kitagawa, T. Best, S. WiU, E. Demler, E. Altman, I. 
Bloch and B. Paredes, Science 327, 1621 (2010). 
M. Greiner, I. Bloch, O. Mandel, T. W. Hansch and T. 
Esslinger, Phys. Rev. Lett. 87, 160405 (2001). 
A. K. Geim and K. S. Novoselov, Nature Mate. 6, 183 
(2007). 

J. P. Reed, B. Uchoa, Y. I. Joe, Y. Can, D. Casa, E. 
Fradkin and P. Abbamonte, Science 330, 805 (2010). 
K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, 
M. I. Katsnelson, I. V. Grigorieva, S. V. Dubonos and A. 
A. Firsov, Nature 438, 197 (2005). 

Z. Y. Meng, T. C. Lang, S. Wessel, F. F. Assaad and A. 
Muramatsu, Nature 464, 847 (2010). 
S. Ladak, D. E. Read, G. K. Perkins, L. F. Cohen and 
W. R. Branford, Nat. Phys. 6, 359 (2010). 
A. H. Castro Neto, F. Guinea, N. M. R. Peres, K. S. 
Novoselov and A. K. Geim, Rev. Mod. Phys. 81, 109 
(2009). 

P. Soltan-Panahi, J. Struck, P. Hauke, A. Bick, W. 
Plenkers, G. Meineke, C. Becker, P. Windpassinger, M. 
Lewenstein and K. Sengstock, Nat. Phys. 7, 434 (2011). 
P. Soltan-Panahi, D. S. Liihmann, J. Struck, P. Wind- 
passinger and K. Sengstock, Nat. Phys. 8, 71 (2012). 
E. Kim and M. H. W. Chan, Nature 427, 225 (2004). 
E. Kim and M. H. W. Chan, Science 305, 1941 (2004). 
A. Cho, Science 336, 661 (2012). 
S. Balibar, Nature 464, 176 (2010). 

E. Kim and M. H. W. Chan, Phys. Rev. Lett. 109, 
155301 (2012). 

[20] G. G. Batrouni and R. T. Scalettar, Phys. Rev. Lett. 84, 
1599 (2000). 



5 



[21] F. Hebert, G. G. Batrouni, R. T. Scalettar, G. Schmid, 
M. Troyer and A. Dorneich, Phys. Rev. Lett. 65, 014513 
(2001). 

[22] P. Sengupta, L. P. Pryadko, F. Met, M. Troyer and G. 

Schmid, Phys. Rev. Lett. 94, 207202 (2005). 
[23] S. Diehl, M. Baranov, A. J. Daley and P. ZoUer, Phys. 

Rev. Lett. 104, 165301 (2010). 
[24] D. Heidarian and A. Paramekanti, Phys. Rev. Lett. 104, 

015301 (2010). 

[25] S. Wessel and M. Troyer, Phys. Rev. Lett. 95, 127205 

(2005) . 

[26] F. Wang, F. PoUmann and A. Vishwanath, Phys. Rev. 
Lett. 102, 017203 (2009). 

[27] R. G. Melko, A. Paramekanti, A. A. Burkov, A. Vish- 
wanath, D. N. Sheng and L. Balents, Phys. Rev. Lett. 
95, 127207 (2005). 

[28] S. Saccani, S. Moroni and M. Boninsegni, Phys. Rev. 
Lett. 108, 175301 (2012). 

[29] L. Seabra and N. Shannon, Phys. Rev. Lett. 104, 237205 
(2010). 

[30] K. K. Ng and T. K. Lee, Phys. Rev. Lett. 97, 127204 

(2006) . 

[31] N. Laflorencie and F. Mila, Phys. Rev. Lett. 99, 027202 

(2007) . 



[32] P. Sengupta and C. D. Batista, Phys. Rev. Lett. 98, 
227201 (2007). 

[33] P. Sengupta and C. D. Batista, Phys. Rev. Lett. 99, 
217205 (2007). 

[34] P. Maher, C. R. Dean, A. F. Young, T. Taniguchi, K. 

Watanabe, K. L. Shepard, J. Hone and P. Kim, Nat. 

Phys. 9, 154 (2013). 
[35] M. Kharitonov, Phys. Rev. Lett. 109, 046803 (2012). 
[36] O. F. Syljuasen and A. W. Sandvik, Phys. Rev. E 66, 

046701 (2002). 

[37] L. M. Duan, E. Demler and M. D. Lukin, Phys. Rev. 

Lett. 91, 090402 (2007). 
[38] S. L. Zhu, B. Wang and L. M. Duan, Phys. Rev. Lett. 

98, 260402 (2007). 
[39] K. Eckert, O. Romero-isart, M. Rodriguez, M. Lewen- 

stein, E. S. Polzik and A. Sanpera, Nat. Phys. 4, 50 

(2008). 

[40] S. Folhng, F. Gerbier, A. Widera, O. Mandel, T. Gericke 

and L Bloch, Nature 434, 481 (2005). 
[41] T. A. Corcovilos, S. K. Baur, J. M. Hitclicock, E. J. 

Mueller and R. G. Hulet, Phys. Rev. A 81, 013415 (2010). 
[42] J. S. Douglas and K. Burnett, Phys. Rev. A 82, 033434 

(2010). 

[43] W. S. Bakr, J. L Gillen, A. Peng, S. Foiling and M. 
Greiner, Nature 464, 72 (2009). 



